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Abstract 

An attractive technique for forming and collecting aggregates of magnetic 
material at a liquid-air interface by an applied magnetic held gradient was 
recently addressed theoretically and experimentally [Soft Matter, (9) 2013, 
8600-8608]: when the magnetic held is weak, the dehection of the liquid-air 
interface has a steady shape, while for sufficiently strong helds, the interface 
destabilizes and forms a jet that extracts magnetic material. Motivated by 
this work, we develop a numerical model for the closely related problem of 
solving two-phase Navier-Stokes equations coupled with the static Maxwell 
equations. We computationally model the forces generated by a magnetic 
held gradient produced by a permanent magnet and so determine the in¬ 
terfacial dehection of a magnetic huid (a pure ferrohuid system) and the 
transition into a jet. We analyze the shape of the liquid-air interface during 
the deformation stage and the critical magnet distance for which the static 
interface transitions into a jet. We draw conclusions on the ability of our nu¬ 
merical model to predict the large interfacial deformation and the consequent 
jetting, free of any htting parameter. 
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1. Introduction 


Synthesis and assembly on the nanoscale is an important goal of con¬ 
temporary science and technology. Magnetic nano/microparticles arise in a 
wide range of industrial and biomedical applications, and so are one target 
for controlled assembly. For exaimle, functionalized magnetic microparti¬ 
cles can be used to separate cells [l| and magnetic microparticles have been 
used in microfluidics for cell sorting, blood cleansing, and magneto-capillary 
self-assembly (see e.g. and references therein). When magnetic nanopar¬ 
ticles such as magnetite are suspended at high concentration in aqueous or 
non-aqueous carrier fluids, the entire system behaves as a continuum of mag¬ 
netic fluid, known also as a ferrofluid. The rheology and interfacial shape of 
ferrofluids can be tuned with external magnetic helds, often in useful ways. 
An example is the application of ferrofluids in adaptive optics that has been 
considered in recent experiments BB- The control of ferrofluid properties 
using magnetic helds also has applications in mechanical sealing and acous¬ 
tics jH , targeted drug delivery B 0, B and treatment of retinal detachment 


Thin liquid dims and droplets are ubiquitous in nature and also appear 
in many technological applications. The understanding of their dynamical 
behavior and their stability is therefore of great importance and has at¬ 
tracted considerable attention in the literature. Recent research into thin 
him and droplet hows has resulted in many experimental and theoretical 
developments, including manipulating him hows via external magnetic or 
electric helds to produce nanoscale patterns. In particular, experiments on 
thin ferrohuid films and drop lets have revealed the formation of a wide range 
of morphologies [13, 11, [ij 13, 14|. Ferrohuids can be manipulated using 


magnetic forces and have been extensively investigated and widely used in 
a variety of engineering applications; see Rosensweig [iB and a more recent 
review by Nguyen [l^. Normal held instability of ferrohuid hlms (and the 
equivalent electric held problem) have been extensively studied in the past, 
see e.g. jl7, 18|. However, despite the increase in the number of applications, 
surprisingly, little can be found in the literature on the direct numerical sim¬ 
ulations of thin ferrohuid hlms in the presence of a nonuniform magnetic held 
(such as is produced by a spherical magnet) and therefore our understanding 
of the instabilities that may occur in these hows is limited. 

An attractive technique for forming and collecting aggregates of magnetic 
material at a liquid-air interface by an applied magnetic held was recently 
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addressed experimentally and theoretically by Tsai et al. [19[. In the exper¬ 
iments described in [l^, a water-based ferrofluid (EMG805, Ferrotec), with 
a density of 1200 kg m“^ and viscosity of 3 mPa s, is suspended in a shal¬ 
low reservoir containing deionized water, with a density of 1000 kg m“^ and 
viscosity of 1 mPa s, to form the magnetic mixture. A spherical permanent 
magnet is then slowly brought close to the magnetic mixture allowing the 
ferrofluid to aggregate and form a static hump at the liquid-air interface (cf. 
Fig. [T]). In these experiments, a distinct boundary that separates the mag¬ 
netic and non-magnetic regions is observed. When the magnet is held suffi¬ 
ciently close to the liquid-air interface, the hump destabilizes and transforms 
to a jet. The theoretical approach developed in [l^ describes a steady-state 
mathematical model for the behavior of the magnetic-particle-laden fluid and 
the particle-free fluid regions. The mathematical model results in jl^ show 
excellent agreement with the experimental data. 

Motivated by this work, here we develop a numerical model for a closely 
related problem: we computationally model the magnetically induced inter¬ 
facial deflection of a magnetic fluid (ferrofluid) and the transition into a jet 
by a magnetic held gradient from a permanent magnet placed above the free 
surface. The system we study differs from that considered by Tsai et al. (l9 


we consider a pure ferrohuid system, while Tsai et al. model a system with 
both magnetically dominated and non-magnetic regions. Figure [T] shows 
a schematic illustrating the setup we consider: the magnetic region occu¬ 
pied by pure ferrohuid, the liquid-air interface, and the spherical permanent 
magnet. The deformation of the ferrohuid-air interface arises as a result of 
the magnetic held gradient induced by the spherical permanent magnet held 
above the huid. We solve the two-phase Navier-Stokes equations coupled 
with the static Maxwell equations. We analyze the shape of the liquid-air 
interface during the deformation stage and the critical magnet distance (from 
the undeformed free surface), for which the static interface transitions into 
a jet. We draw conclusions regarding the ability of our numerical model to 
predict the large interfacial deformation and the consequent jetting, free of 
any htting parameters. The numerical model provides a realistic and accu¬ 
rate framework for predicting the evolution of magnetic liquids based on the 
Navier-Stokes equations. 

We describe the details of the numerical model in Sec. |2j In Sec. |3l 
we describe a numerical boundary condition that may be implemented to 
simulate non-uniform magnetic helds. In Sec. HI we present the numerical 
results and the comparison with experimental observations. In Sec. 0 we 
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Figure 1: A schematic illustrating the computational setup and the coordi¬ 
nate system used. A spherical permanent magnet of radius Rm is centered 
at distance L from the initially undeformed film (red dashed line), which has 
a depth Lq. The magnetic force deforms the interface, Shi/, into a hump 
(red solid line). A unit normal outwards from the interface is denoted by n. 
A typical computational domain, ff, and its boundary, dVt^ is shown by the 
dash-dotted line. Due to symmetry, only half of the domain is simulated. 


give an overview and future outlook for improving our modeling. 

2. Mathematical Model 

Here we briefly describe the theoretical models that serve as a basis for the 
proposed numerical studies. The coupled motion of a ferrofluid surrounded 
by a non-magnetic fluid is governed by the (static) Maxwell equations, the 
Navier-Stokes equations, and a constitutive relationship for the magnetic 
induction B, magnetic field H, and the magnetization M. The magnetostatic 
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Maxwell equations for a non-conducting ferrofluid are, in SI units, 


V-B = 0, VxH = 0, B(x, t) 


/x/H in ferrofluid 
/XmH in matrix. 


where fif denotes the magnetic permeability of the ferrofluid and is the 
permeability of the matrix fluid. For our application, the matrix fluid is air, 
which has a permeability very close to that for a vacuum, /Iq. Therefore, we 
shall consider pm = po throughout this mathematical model. A magnetic 
scalar potential tfj is dehned by H = and satishes 


V ■ (/iVV') = 0, 


( 1 ) 


where p = /io and p/ in the matrix and ferrofluid, respectively. We will 
assume that the magnetization is a linear function of the magnetic held given 
by M = XH, where y = (/i///io — 1) is the magnetic susceptibility j^. The 
magnetic induction B is therefore B = /ro(H -l- M) — /io(l + y)H. 

The huid equation of motion is described by the conservation of mass and 
momentum (Navier-Stokes) equations 


P 


V-u 



+ (u-V) 



0 , ( 2 ) 

-Vp -h V ■ (2pD) -h V • -t- Fs -1- pg, (3) 


where F^ denotes the surface tension force per unit volume (presented as a 
body force |2]|), p is pressure, u is velocity, D = |(Vu-|- (Vu)^) is the rate 
of deformation tensor (where T denotes the transpose), p is viscosity, p is 
density, Tm is the magnetic stress tensor, and g is the gravitational accelera¬ 
tion. The total stress is r = —pi -|- 2pD -|- r^, where I denotes the identity 
operator. The magnetic stress tensor of an incompressible, isothermal, mag¬ 
netizable medium is 22 


T^ = -^Hn + pmi^ 


where H = |H|. These equations must be solved subject to suitable boundary 
and initial conditions, discussed in Sec. |3] below. 


3. Numerical Methodologies 

We will use an Eulerian framework, where the material moves through a 
stationary mesh, and therefore a special procedure will be required to track 
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the interface between fluids. We will use the volume of fluid (VOF) method 


to track the interface between the ferrofluid and the matrix fluid 


24, m 


In this way, the VOF formulation describes each fluid by assigning a volume 
fraction function, f{r,z,t), as 


f{r,z,t) = 


1 in ferrofluid 
0 in matrix 


(4) 


(see Fig. [T]for the polar coordinates {r,z) used here). The position of the 
interface is computed from f{r,z,t) by reconstructing the curve where the 
step discontinuity takes place. In this work, the reconstruction is a ‘piecewise 
linear interface calculation’ (FLIC) j^. To track the interface, f{r,z,t) is 
advected by the flow 

^ + V.(u/) = 0. (5) 

In the Navier-Stokes equations (Eq. ([3])), includes the continuum body 
force due to interfacial tension 


Fs = cTKh^n, (6) 

where the unit normal n = V//|V/| and 6s = |V/| is the delta function at 
the interface. 


3.1. An Optimal Method for Simulating Non-Uniform Magnetic Fields 
Special care must be taken when computing the solution of the Maxwell 
equations to account, accurately and robustly, for the non-uniformity of 
the magnetic held. The boundary condition on the magnetic held is re¬ 
constructed from an existing solution of the magnetic held due to a spherical 
magnet, in the absence of the ferrohuid. We let ip a denote the magnetic po¬ 
tential for this inhnite-domain analytical solution, satisfying Laplace’s equa¬ 
tion, in cylindrical coordinates with azimuthal symmetry, 

1A (r^\ + ^ = 0 

r dr \ dr J dz"^ 


We assume a magnetic potential generated by a uniformly magnetized spher¬ 
ical magnet, as in 1^, hence 


'^a{r,z) 


MrnRljL - {z - Lo))r 

3 |(L -(z- + N'" 


+ constant. 


(7) 
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where is the magnetization of the magnet and Rm is the radius of the 
magnet. The magnet center is at (r, z) = (0, L + Lq), where L is the magnet 
distance from the undeformed interface and Lq is the initial unperturbed in¬ 
terfacial thickness in our problem setup (see Fig. [1]). The analytical magnetic 
field Ha = is thus 


M^KiQz - Lq) - L)r ) f [2((z - Lq) - Lf - r^] 1 , 

[(L - (2 - L„))= + rf/d \ 3 [(L - (2 - L„))= + 


where r and z denote the unit vectors in the r and z directions, respectively. 
This analytical solution, Eq. ([7]), provides the boundary conditions that we 
will impose on the magnetic potential at the boundaries of our computational 
domain. 

The numerical magnetic potential 'ip is calculated by solving Eq. o on 
the computational domain hi (see Fig. [T]) employing the boundary conditions 
obtained using the analytical solution 'ipa{,r,z) above. In all simulations pre¬ 
sented, the spherical magnet lies outside the computational domain. The 
numerical solution to the elliptic partial differential equation (Eq. ([1])) will 
be complicated by the fact that the magnetic permeability, /i, experiences a 
jurnp across the interface [^. We will use a multigrid algorithm, as described 
in |20|, to solve Eq. ([T]) subjected to the boundary conditions described below. 

The boundary conditions for ip on the domain boundaries dQ are defined 
as 


dip 

dn 


dipa 

dn 


on dQ, 


(9) 


where d/dn = nr ■ V, and np denotes the normal to the boundary dQ. In 
order to impose the boundary condition in our numerical model, we perform 
a transformation of variables to C: ip = ipa + C, where ipa is the potential field 
without the magnetic medium, Eq. ([71). One can then rewrite Eq. ([T]) such 
that 

V-(/iVC) = -V-(/iVV'a), (10) 

where V ■ (/iV^’a) vanishes everywhere except on the surface between the 
ferrofluid/matrix interface dVtf and 

^ = 0 on dVt. 
dn 
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( 11 ) 





Though we do not use this in the following analysis, we note that our 
numerical approach may also be used to model the nonlinear magnetic sus¬ 
ceptibility of ferrofluids; for example, the well-known Langevin function 
L{a)=coih.a — a~^ can be implemented to describe the magnetization be¬ 
havior of the ferrofluid versus the strength of the magnetic held H, 

M,„, . (MS) H , 

where the saturation magnetization Mg and the magnetic moment of the 
particle enter as parameters, T denotes the absolute temperature, and ks is 
the Boltzmann constant. 

3.2. Pressure, Velocity, and Volume Fraction Boundary Conditions 

We impose symmetry boundary conditions at r = 0. At the top, bot¬ 
tom and right boundaries, we impose solid wall boundary conditions for the 
pressure and the velocities. The boundary condition for the volume fraction 
function at the top and bottom walls is that / = 0 and for the right wall is 
that the interface has zero slope ^ = 0. 


4. Results and Discussion 


Numerical simulations are presented in three parts. First we present tests 
of the numerical implementation by focusing on the resulting magnetic helds 
and by studying the convergence of the solution for the magnetic held with 
grid rehnement. Second, we present several numerical solutions of our system 
that demonstrate both the steady dehection regime, and the transition from 
this regime to the unstable jetting regime. A dimensional analysis in fl^ 
shows that the steady interfacial deformation can be characterized by the 
magnetic Bond number that captures the dominant balance between the 
magnetic and the surface tension forces. Following this approach, we present 
a parametric study by varying the magnetic Bond number within the steady 
interface dehection regime. 

Our numerical setup closely follows that in 
magnet of radius R 
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a spherical permanent 
and magnetization = 10® A m“^ is placed at a 
distance L from the initially undeformed ferrohuid him (see Fig. [1]). The fer- 
rohuid has a density Pf = 1200 kg m“^, viscosity pf = 3 mPa s, and surface 
tension a = 0.07 N m“^. The ferrohuid is assumed to have a constant mag¬ 
netic susceptibility, y, which will be determined by comparing the numerical 
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Figure 2: The computational domain, hi, and the computed magnetic held 
strength, H, with a drawing of a permanent magnet overlaid. The spherical 
magnet lies outside the computational domain. Number of grid points in the 
2 ;-direction is Nz = 240 and in the r-direction is Nr = 120. 


result with the experimental data in 19| (note that this value is often not 
reported by manufacturers so is determined by users). 


4.1. Magnetic Field and Imposed Numerical Boundary Condition 

We hrst demonstrate the effectiveness of our methodology for computing 
the magnetic held by presenting the results of the simulated applied magnetic 
held using the numerical boundary conditions described in Sec. 13.11 The 
following relative error norms are dehned 


11“ I l 2 ,A] 2 iax 

where 


1771 


oo.A 


1771 


oo,A 


max 


1771 


oc-^max 


1771 


2,A 




and 1177| |oo,A = Max (|77jj|) 


A ’ 
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Number of grid points in z-direction, 


Figure 3: The convergence of the L 2 and L^o relative error norms between 
the H{r,z) computed at A = Lz/Nz{= Lr/Nr) and at Amax = -f^2/240, for 
the computational setup illustrated in Fig. [2l 


where i and j are indices of a computational cell. Figure |2] shows the compu¬ 
tational domain r G (0, L^) and G (0,L^), where {Lr,Lz) = (1.5 mm, 3 
mm), with a 1 mm radius permanent magnet placed 5 mm above the r-axis. 
The mesh size A = Lr/Nr = L^/N^, where Nr and N^ are the number of 
grid points in r and z directions, respectively. The computed magnetic held 
strength, H, is also shown in Fig. [2j 

Figure |3] exhibits the convergence of the numerical method with spatial 
resolution for computing the magnetic held. As illustrated, a second-order 
convergence is obtained for the mesh rehnement for both L 2 and L^o error 
measures. The hgure also reveals the smallness of the errors even at a coarse 
grid resolution. Finally the accuracy of the numerical solution is also as¬ 
sessed by comparing the computed magnetic held with the exact solution for 
a spherical magnet in an unbounded region given by Eq. (|8]); the maximum 
errors in H are of the order 0.01%. These results conhrm the overall con¬ 
sistency and accuracy of the computational scheme with the implemented 
numerical boundary conditions as described in Sec. 13.11 


4 . 2 . Interfacial Deflection and Transition to Jetting 

Tsai et al. describe the results of experiments and mathematical analysis 
of the deformation of a free surface by a magnetic force from a spherical per¬ 
manent magnet 1^. As mentioned above, in that work the system consists 
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of both magnetic and non-magnetic regions such that the magnetic particles 
collect at the interface, causing deformation of the free surface to form a 
hump; and when the magnet is brought sufficiently close, the hump transi¬ 
tions to a jet. Within the non-magnetic region, the pressure drop across the 
static interface, and hydrostatic stresses, are accounted for by the Young- 
Laplace equation; within the magnetic region, the liquid is treated as a static 
ferrofluid that has a homogeneous magnetic susceptibility, and the modihed 
form of the Young-Laplace equation (modihed to account for the huid mag¬ 
netization) describes the pressure drop across the static interface. The steady 
state of the interface shape is then determined by two second-order differ¬ 
ential equations solved by imposing appropriate boundary conditions. By 
exploiting the smallness of the vertical dehection of the huid surface com¬ 
pared to the horizontal extent occupied by the paramagnetic particles, and 
the smallness of the magnetic susceptibility, Tsai et al. also analyze their 
mathematical model asymptotically in the non-jetting regime. Using the 
asymptotic solution and a scaling argument, they compute the maximum in¬ 
terfacial dehection at the center of the hump to predict when the transition 
to jetting occurs 
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The system we investigate computationally here dihers from that of [19 


because our liquid region contains a homogeneous ferrohuid. Despite this dif¬ 
ference, here we also observe the same qualitative features: when the magnet 
is positioned sufficiently far from the ferrohuid, the deformation of the in¬ 
terface is static; and when we reduce the distance between the magnet and 
the undeformed interface, the deformation destabilizes into jet formation. In 
Fig. m we present the results of direct computations of our system in the 
static deformation regime. In these simulations, the computational domain 
is Lr = 15 mm by = 2.5 mm and the undeformed him is initially located 
at Lq = 1 mm. A magnet of radius Rm = 3.2 mm is placed at distance 
L = 6.33 mm (Fig. 0^) and 5.65 mm (Fig. |3 )d) from the undisturbed him. 
In Fig. m we also reproduce the experimental results of [l^ to compare and 
contrast the two systems. 
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(a) 



r (mm) 

(b) 

Figure 4: Experimental data (not pure ferrofluid) from jl^ shown alongside 
our computed steady-state ferrofluid-air interfacial profiles for (a) L = 6.33 
mm and (b) L = 5.65 mm. Black solid lines have y = 7.5 x 10“^, chosen 
to match maximum free surface height to data; green solid lines have y = 
5 X 10“^, chosen to £t the data as a whole. The results illustrate the height 
of the hump and the radial spread of the ferrofluid-air interface. 

We use the comparison of the computed steady-state hump height, h, with 
the experimental measurement to determine the magnetic susceptibility of 
the ferrofluid. This comparison gives y = 7.5 x 10“^, a value consistent with 
the hnding in [l^. In Fig. 01 we also demonstrate the results when using 
y = 5 X 10“^, chosen to £t the data as a whole rather than just htting to 
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L = 6.33 mm 
L = 5.65 mm 



Figure 5: Static ferrofluid-air interfaces when varying L; for L < 5.05 mm, 
no static hump can form. A decrease in L results in decreasing the radial 
spreading. The hump also becomes less pronounced with increasing L. 


the maximum hump height. Although the results demonstrate consistency 
between the direct numerical solutions and the experiments in capturing the 
height of the hump, the agreement in the prohle is not quantitative, owing 
to fact that the experimental system is not pure ferrofluid. 

With this caveat in mind, we next investigate the transition to jetting 
with y = 7.5 x 10“"^. In Fig. |5l we show a series of static ferrofluid-air 
interface prohles as the distance L between magnet center and undeformed 
ferrofluid interface is reduced. As the magnet separation L is decreased 
incrementally, we hnd that there exists a threshold for the jetting transition: 
for L < 5.05 mm, we no longer see the free surface evolving to a static 
conhguration; instead, the ferrofluid-air interface destabilizes and forms a 
jet. The maximum sustainable steady interfacial deflection that we were able 
to compute occurs at L = 5.05 mm, consistent with the experimental hnding 
in [l^. Also consistent with the experimental observations, we hnd that 
the steady interfacial dehection increases and the radial spreading decreases, 
with decreasing L, until the hump transitions to jetting. In Fig. |6l we plot 
interface evolution when the magnet is sufficiently far away {L = 5.25 mm) 
that the interface reaches a static state (Fig. [6](a)); and when it is close 
enough (L = 4.95 mm) that jetting occurs (Fig. [Hl^b)), where no sustainable 
static interfacial dehection can be obtained. Interestingly, for L = 5.25 mm 
(Fig. in](a)), we hnd an overshoot in deformation; i.e., an initial elongation 
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Figure 6: The ferrofluid-air interface, (a) t = 0,25,50,100,200,450 ms 
(from bottom to top) for L = 5.25 mm. The interface deforms and reaches a 
static state (green dashed line); the interface undergoes a transient overshoot 
before settling down to a static shape, (b) t = 0, 5, 25, 50,100,150, 200, 248 
ms (from bottom to top) for L = 4.95 mm. No sustainable steady interfacial 
deflection can be obtained and the dynamic interface grows until it touches 
the top boundary of the computational domain at 250 ms (not shown here). 
The magnet radius Rm = 3.2 mm and the experimentally htted value of the 
magnetic susceptibility y = 7.5 x 10“^. 

followed by a retraction to a static state. We attribute the overshoot to 
the competition between magnetic, capillary, and inertial effects. Also, as 
illustrated, for L = 4.95 mm (Fig. EJ^b)), the interface becomes unstable 
and stretches until it touches the top boundary of the computational domain 
while for L = 5.25 mm, a stable interface is achieved. We also note that, 
when the interface destabilizes, the deflection increases very rapidly as the 
ferrofluid-air interface approaches the magnet. 

We also investigate the effect of increasing the magnetic susceptibility and 
the magnet size. In Fig. [71(a), we plot the ferrofluid-air interface for y = 1.75 
as a function of time for L = 6.75 mm and Rm = 1-5 mm. We also plot the 
case when y = 0.75 in Fig. [71(b). These results indicate that increasing the 
magnetic susceptibility can destabilize a previously stable interface, moving 
it into the jetting regime; and that when the magnet size is reduced, a much 
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(a) (b) 

Figure 7: The ferrofluid-air interface, (a) t = 0,12.5,25,37.5 ms (from 
bottom to top) for x = 0.75. The interface deforms and reaches a static state 
(green dashed line) at around t = 150 ms. (b) t = 0, 2.5, 5, 7.5,10,12.5,15 ms 
(from bottom to top) for x = 1.75. The interface touches the top boundary 
of the computational domain at around 16 ms (not shown here). L = 6.75 
mm and Rm = 1.5 mm. The flat interface shows the initial position of the 
ferrofluid-air interface. 


higher magnetic susceptibility is required to destabilize the interface into 
jetting. These predictions can be used to hnd the optimal magnet size for 
destabilizing the interface based on the susceptibility of the magnetic fluid. 
For a hxed magnet distance L and magnet size Rm, transition to jetting 
occurs for a certain value of the magnetic susceptibility x- When the magnet 
is brought slightly closer (decrease L), the transition occurs for a smaller x 
value. Figure [H] shows the results for increasing the magnetic susceptibility, 
which leads to an increase in the interfacial deflection for L = 6 mm and 
Rm = 1-5 mm. This hgure shows that, for a sufficiently high magnetic 
susceptibility, the surface tension of the ferrofluid can no longer sustain the 
deformation and jetting again occurs. We find that in this case, for x ^ 0.1, 
no static hump can form. We also note that for x = 0.75 in Fig. [HI a secondary 
instability appears to form at the tip of the jet, although this instability will 
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Figure 8: The ferrofluid-air interface when varying y; for x > 0.1, no static 
hump can form. Increasing y results in transition to jetting. For y = 0.75, 
the dynamic interface stretches until it touches the top boundary of the 
computational domain, while the interface reaches a static state for y = 0.075 
and X = 0.1. In all cases L = 6 mm and Rm = 1-5 mm. A secondary 
instability (or necking that will grow into a pinch-off structure) appears to 
develop at the tip of the jet for x = 0.75. 

not have sufficient time to grow (most likely into a pinch-off structure) before 
the tip touches the magnet. 

4-3. Scaling Model for Interfacial Deflection 

Finally, we investigate the maximum interfacial deformation, h(0), de¬ 
fined as the maximum height of the hump at r = 0, as L varies. To char¬ 
acterize this behavior, it is helpful to represent the results in terms of the 
following dimensionless parameters [l^: the magnetic Bond number 

_ Xt^oM-rn^rn 

— 1 
a 

which represents the ratio of the magnetic to surface tension force; and 

L* = and h* = 
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Figure 9: The maximum interfacial deformation, h*{0), versus the non- 
dimensional distance between the magnet center and the initially flat him, 
L*, for various magnetic Bond numbers, Bom- Rm = 3.2 mm and a = 0.07 
N m“^. Critical points are marked by red arrows; a smaller L* will result in 
no static prohles. 


which give the dimensionless distance between the magnet center and the ini¬ 
tially undeformed ferrohuid him; and the dimensionless height of the hump, 
respectively. In Fig. |9l we present the results of the maximum dimension¬ 
less interfacial deformation, h*(0), predicted by our computational model, 
versus dimensionless magnet distance L*, for several values of Bom,. Again, 
we see that as the magnet is successively brought closer to the interface, the 
maximum hump height increases. We also present the ehect of the magnetic 
Bond number, showing an increase in the maximum interfacial deformation 
as the magnetic Bond number is increased. For all values of Bom considered, 
the results show that the maximum interfacial dehection does not vary sig- 
nihcantly prior to jetting, while the magnitude of the maximum deflection 
increases rapidly for smaller L*. Past the critical values of L*, no static pro- 
hle can be obtained (the critical points are shown by red arrows in Fig. IH]). 
In general, the results indicate that when using a ferrohuid with a smaller 
susceptibility, the magnet must be placed closer to the interface to produce 
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Figure 10: Non-dimensional maximum interfacial deformation, h*{0), plotted 
against for the results in Fig. |9l The solid line is the scaling 

£t, h*{0) oc (Bo™)1/VL*T 


jetting. 

Analogous to the study in 
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we identify the natural scaling for the 
characteristic interfacial deflection he by 


^2 _ hoXMlRl, _ BOmRl, 

representing the dominant balance between the magnetic and surface tension 
forces. Numerical results in Fig. |9] suggest that, prior to jetting, the maximum 
interfacial deflections do not vary signihcantly for different values of Bo^- 
Motivated by this observation, we use Eq. da to determine the scaling 
relationship between the maximum interfacial deflection and the parameters 
that we vary in our numerical simulations: L* and Bom- In Fig. [101 we plot 
h*(0) versus \/Bo^/L*^, predicted by the scaling model Eq. flT3|) . Further 
from the transition, the simulations results collapse remarkably well onto 
this single scaling relationship. Close to the jetting, however, we observe a 
deviation from the scaling law, whereby the assumption of slowly varying 
h*(0) versus L* breaks down. 
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5. Conclusions 

In this paper, we have carried out direct computations of the deformation 
of a ferrofluid-air interface under an external magnetic field gradient gener¬ 
ated by a spherical magnet placed at a fixed distance from the interface. 
We showed that, for sufficiently large distances, the free surface deforma¬ 
tion is static, while below a threshold value, the interface destabilizes and 
forms a jetting fluid. These features are also observed experimentally, in a 
somewhat different system in which the ferrofluid is suspended in water (l^ . 
Our numerical results provide data that can be used to determine the max¬ 
imum deflection of a ferrofluid in the presence of a magnetic held gradient 
generated by an external permanent magnet and the consequent transition 
to jetting. Additionally, a simple scaling law allowed us to collapse our nu¬ 
merical results from a series of conhgurations onto a single power law. It is 
interesting to note that, while the key features observed in the experiments 
are reproduced in our pure ferrohuid simulations (i.e., static interfacial de¬ 
formations at large magnet separation or low magnetic susceptibility, with 
transition to jetting at small magnet separations or high magnetic suscepti¬ 
bility), the static interfacial prohles obtained in the two systems are rather 
different. In particular, the interfacial prohles for the water-ferrohuid mix 
are more “peaked” under the same conditions than those for pure ferrohuid; 
see Fig. |U This could be due to the water providing a lubricating “slip” layer 
that allows the ferrohuid to be more mobile in the experiments. Investigating 
such mixed systems computationally will be deferred to a future publication. 
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